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Abstract. We develop a theoretical analysis of four- wave mixing used to generate 
photon pairs useful for quantum information processing. The analysis applies to a 
single mode microstructured fibre pumped by an ultra-short coherent pulse in the 
normal dispersion region. Given the values of the optical propagation constant inside 
the fibre, we can estimate the created number of photon pairs per pulse, their central 
wavelength and their respective bandwidth. We use the experimental results from a 
picosecond source of correlated photon pairs using a micro-structured fibre to validate 
the model. The fibre is pumped in the normal dispersion regime at 708 nm and 
phase matching is satisfied for widely spaced parametric wavelengths of 586 nm and 
894 nm. We measure the number of photons per pulse using a loss-independent 
coincidence scheme and compare the results with the theoretical expectation. We 
show a good agreement between the theoretical expectations and the experimental 
results for various fibre lengths and pump powers. 
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1. Introduction 

Quantum information can be encoded in various media including atoms, ions and 
electrons. It is photons however that are the most useful for transporting quantum 
information between separate locations. They are often called flying Q-bits and they 
appear to be the fundamental resource for quantum communications experiments PHEJ- 
Their robustness to decoherence also leads them to be involved in various multiphoton 
and linear optical logic experiments [31 HJ Ej. One of the main resources needed for 
these quantum optics experiments is the ability to create these photons in pairs. When 
simply correlated in time they can be used to create heralded single photons |SHZj (closely 
approximating true single photon sources). More interesting is when these photons have 
strongly correlated properties; they are entangled [HIE]- Such sources continue to be 
exploited in fundamental multiphoton experiments such as quantum teleportation [TOJ 
or to create cluster states [TTJ [T2J H3] • 

The preferred sources for such experiments until recently have been three-wave 
mixing in x non-linear birefringent crystals (HI E|. These sources are inherently 
wide band, low brightness (per nanometer, per single mode) sources. The increasing 
number of photons involved in current experiments is requiring brighter sources and 
narrower photon bandwidths. Hence, periodically poled fibres [TJ] and periodically 
poled waveguides of lithium niobate ^Hj have been shown to be bright pair photon 
sources. In poled fibres ^H] the low non-linearity limits the brightness while in planar 
waveguides the non-circular mode limits the coupling efficiency into optical fibres |17j . 

It is well known that parametric gain can arise from the x^ non-linearity in optical 
fibres [T%1 IT!?] and phase matching can be achieved by using the modulation instability 
when pumping fibres in their anomalous dispersion regime. Various pair photon 
generation experiments have been performed in this regime [201 EU 1221 1231 121] • The 
photon pairs are generated close to the pump wavelength and are always accompanied 
by a significant Raman background and careful filtering is required. We have shown that 
phase matching can be obtained for widely spaced wavelengths by pumping photonic 
crystal fibre (PCF) close to the zero dispersion wavelength in the normal dispersion 
regime [23 [23 123 EH] • Recently, we have used a fibre with the zero dispersion point 
in the near infra-red (715 nm) and pumped with a picosecond pulsed laser, red detuned 
a few nanometers into the normal dispersion regime, to demonstrate a high brightness 
source of photons pairs [2Hj . 

We are seeking to develop a source which may be applicable for future quantum 
interference experiments involving three or more photons created as two or more pairs. 
Interference effects between separate pairs of photons can be studied by overlapping 
photons with a time uncertainty shorter than their inverse bandwidth or coherence 
length [30] • This restricts us to sources pumped by ultra-short laser pulses where the 
bandwidth is of order nanometers and also requires a high efficiency of collection. 
Theoretical models of correlated photon pair generation via four-wave mixing have 
already been investigated in the CW approximation JH], but this approximation no 
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longer stands if we want these photons ready for quantum interference experiments. We 
propose here the first quantum model to describe four-wave mixing to generate photon 
pairs in the pulsed regime in a single mode fibre. 

In this paper, we investigate the generation of correlated single-photon pairs in 
photonic crystal fibres from both theoretical and experimental points of view. We first 
describe the four-wave mixing process in the quantum mechanical regime. We then 
analyze and compare the experimental results from correlated photon pairs generated 
in a photonic crystal fibre to the aforementioned theory. The paper is structured along 
the following lines: sectional presents the quantum theory of four- wave mixing in optical 
fibres, while section El is dedicated at tagging experimental parameters to the theory, 
section 0] describes and analyzes the experimental results from our photonic crystal 
fibre. We investigate in this section the agreement between theoretical expectations and 
experimental results, section El outlines the possibilities of using this fibre as a source 
of multi-photon pairs for quantum information processing. Finally, we briefly conclude 
in section |B1 



2. Four-wave mixing process in fibre with pulsed pumping 

2.1. Definition of the involved fields 

We are looking at the process of four wave mixing (FWM) based on the non-linear 
susceptibility x m a photonic crystal fibre where two pump-photons are converted 
into a pair of signal and idler photons. We are specializing to energy matched situations 
where 2Q p = uj s + cjj. The related interaction Hamiltonian in a volume within the fibre 
is: 

H mt = J v UdV (1) 

and the energy density associated with four fields coupled in a third order non-linear 
medium is: 

U = e x (3) (fl p> uj s , u i )E%E.E i (2) 

Here eo is the permittivity of free space allowing us to measure x m units of V~ 2 m 2 . 
In the fibre optic community, it is common to rather use the nonlinear refractive index 
n 2 defined as follow: 

where uq is the refractive index of the fibre. 

For a nonstationary pumping pulse, the field amplitude can be decomposed into its 
Fourier components in the following form: 

/+oo 
£ p (uj,x,y,z,t)duj (4) 
-oo 

In our guided configuration all the involved fields are propagating in the same 
direction, hence we insert e p (r) as the transverse spatial variation with normalisation 
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/ l e p(^*)| 2 drdO = 1 inside all our equation and carry out the calculation along the Z-axis 
which is the fibre axis. We assume that E p (x, y, z, t) is a strong classical gaussian pump 
beam of full- width half maximum bandwidth Au p and central frequency Q p . We focus 
on one of its Fourier component defined as follow: 

£ p (cu) = ^G P (cu) (e-iWp+")t-(k P -vP P )z) + e ;[(n p +^)Mfc p - 7 p p ) 2 ]^ ^ ^ 

where •yPp is the self-phase modulation coefficient acquired by the pump pulse as it 
propagates within the fibre; G p (u) = e ~ is the weight of the component Q p + uj of 
the pulse and the variance characterizing the pump pulse. The variance is linked to 

the bandwidth of the pulse by Aj p = 2 v^ n ( 4 ) . We define the energy density within a 
pulse so that (Parseval's theorem): 

£ p (u)S;(u)du = e / -^e'^ du (6) 

While the strong pump pulse remains classical, here the two-photon modes associated 
with E s and Ei are quantized using the creation and annihilation operators and we 
use similar spatial terms e/(r) to describe their transverse variation with normalisation 
J \ei(r)\ 2 drd6 = 1. 



Ei = E (a\e-*" - a t e^) e^r) (7) 

where £ 3 defines the quantization volume. 
2.2. Interaction Hamiltonian 

We calculate the interaction Hamiltonian for one monochromatic component of the 
pump pulse and will sum up over all the possible combinations between the components 
at a later stage. Starting with the interaction Hamiltonian for the parametric process 
and following the steps in Appendix A , we can easily derive the interaction Hamiltonian 



for two pump-components u p i and uj P 2 as: 



1 



'Hint — '"'^V a2 (^P 1 ) 



CO 



p2) 





+ h.c. 









^ila]e- i{2np+ulpl+u)p2)t Lsinc 

where I = jj e p i(r^e P 2{r)e s (r)ei(r)drd9 is a normalised factor; Ak = k p i + k p2 — k s — ki- 
2'yP p is the phase matching coefficient and S is the gain parameter defined as follow: 



(8) 




The sine function in (jSJ) gives rise to the phase-matching condition between the pump 
components, signal and idler fields. The energy associated with the individual signal 
and idler fields in the fibre is defined by 

H = h[u s (a\a s + 1/2) + 00^% + 1/2)] (10) 
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leading to a total Hamiltonian to be: 

H to t = H o + Hint (11) 

Remembering that this Hamiltonian applies only for two arbitrary monochromatic pump 
components, we will later integrate over the full bandwidth of the pump pulse. 



2.3. The Heinsenberg equation of motion 

In the standard Heisenberg representation the evolution of any operator A is defined by 
^ = - n [A,H M \ (12) 
In the case of A = a\ and assuming ai(t) = di(t)e~ luJlt , equation gives: 

al(t -> oo) = 4(0) ^G p (uj pl )G p (u P 2)ai(0)5(Au)sinc(—!^J (13) 

When integrating to obtain (fT3*j) we assumed that the system (nonlinear medium and 
the pump components) is turned on at time zero and runs continuously. This gave 
rise to the ^-function in Au for the frequencies of the fields to guarantee the energy 
conservation at single-photon level for a monochromatic pump component. 

We then perform the integration over all the possible combinations of ou p i,u P 2 
within the pump bandwidth. The full details of the following calculation are given 
Appendix B but basically if we define the interaction satisfying the basic energy 



in 



conservation 2Vt p = lu so + uji and introduce AuJk s as the frequency shift for the signal 
mode k s compared to the central one u SQ (respectively Acu fci for the idler mode ki), the 
5-function then helps to carry out the double integral over {uj p i,lj P 2} and for aj(oo) we 
get 



aj(oo) — "'i 



SIL v^F „ |a^ s+ a^]V _ (AkL\ 



^ (0) ~*27Z a ? C >^M0)«nc[—) ( 14 ) 

It is interesting to note that in (fPfj) the dirac function (standing for the energy 
conservation in monochromatic pump regime) turns here into an exponential function 
guaranteeing the energy conservation in the broadband pump regime. 

In order to calculate the mean number of photons (N s ) over all the possible modes 
k s , we define the operator W s (t) = -7= J2k s so that 

(N.) = ,{0\,(0\W*W,\0)Mi 

= 2 h£\ 2 JJ sinc > e - ^ j^SL ( 15 ) 

where v 9s and v g . are the group velocity of signal and idler wave in the fibre. 

We show in the following section that (|15|) becomes easy to integrate under certain 
approximations and we derive the expected mean number of photon per pulse at the 
output of our photonic crystal fibre. 
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3. From theory to experiment 

Figure H shows the microstructured fibre used in our experiment. The fibre has a core 
diameter of 2 fjm and presents a zero dispersion wavelength A = 715 nm. Typical 




Figure 1. Electron microscope image of the PCF used with core diameter d w 2 /jm 
and Aq = 715 nm. 



spectra obtained when pumping with 708.4 nm light and aligning the polarization on 
one of the axes of the fibre are shown in figure El The asymmetry of the micro-structure 
makes the fibre slightly birefringent. As a consequence, when we do not launch into 
the correct fibre polarization axis, we see two peaks corresponding to different phase 
matching conditions for the different fibre axes. The absence of the second peak here 
implies our source is strongly polarised. 

The short wavelength sits at 586.4 nm and the corresponding idler at 893.9 nm. The 
spectra were taken at a pump power of 1 mW and flat background comes entirely from 
the electronic bias errors in the CCD. Then by tuning the pump photon wavelength and 
monitoring similar spectra, we are able to observe the evolution of the signal and idler 
wavelength versus the pump wavelength (see figure EJ). These wavelengths are those 
which satisfy the following energy conservation: 



2f2p — ui s 



UJi 







and phase matching equations 

. n„ f2„ n s lu s rii Ui 



2 7 P = 



(16) 



;i7) 



where n p s i are the refractive index of the medium at pump, signal and idler wavelengths. 
The phase shift acquired by the pump waves (self-phase modulation) is defined by the 
factor 7P [IB]. We also present on figure El the theoretical phase matched curve using 
the propagation constant for a simple strand of silica in air as a close approximation [3T] 
to the fibre used in this experiment. This approximation does not give a good agreement 
between the theoretical expectations and the experimental results as it does not take 
into account the role played by the honeycomb structure. 



3.1. Phase matched solutions 

As the percentage of silica in the cladding is not negligible, the accuracy of this model 
is not perfect. To solve this problem we use a novel method [22] to determine the 
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Figure 2. Fluorescence spectrum of the signal, pump and idler photons. The 
measurement was integrated over ten seconds. The number of counts is proportional 
to the number of photons detected by the cooled camera. The photon spectra are 
respectively centered at 586.4 nm, 708.4 nm and 893.9 nm. They feature an associated 
FWHM bandwidth of 4.5 nm, 0.4 nm and 9.6 nm. The small peak at 598 nm is 
attributed to unwanted background light or a detector fault since it is visible even 
when the laser is blocked 



propagation constant of a photonic crystal fibre. This method models the PCF as 
a step-index fiber with a core index set to that of fused silica, and a cladding index 
set to the mean index of the PCF's cladding. We used an effective cladding index of 
n = 1.05 (corresponding to 90% air-filling fraction) to retrieve the propagation constant 
of the fibre and obtain the fit shown in figure El There is a very good agreement, 
confirming that our core diameter is actually 2 fjm and that the honeycomb structure 
does contribute to the effective refractive index of the fibre cladding, perturbing the 
propagation constant from that generated by the naked strand model. 

We have now a good knowledge of the propagation constant inside our fibre. This 
numerical data will allow us to simplify the calculation of (|15jl and therefore predict the 
bandwidth and the mean number of photons per pulse. 

3.2. Number of photons per pulse and associated bandwidth 

The integral from (fTK|) can be sorted out if we now develop the Afc-function around the 
phase matched frequencies (2k po — k so — k io — 2jP=0) using a first order Taylor power 
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Figure 3. Nonlinear phase-matching diagram for the process 2fl p — > u> s + Wj. The 
dashed curve is the simulation using a plain strand of silica suspended in the air, while 
the straight curve is the fit used to retrieve the actual refractive index of the fibre. 
The black dots stand for the experimental data while their error bar is proportional 
to their measured FWHM. 



series. 



Ak 



2k 
2 
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dk 
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Alu p — 



2 7 P 
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duj 
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dkj 
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With the help of energy conservation, we can remove the Au p variable and find: 
Ak = (M s - M p ) Auo ks + {Mi - M p ) Auo ki 



(19) 



(20) 



where Mi 



1 dw 



+ ni 



We now have the products of two functions depending only on Au; fcs and Auj k . : The 
first one is smc 2 ^^|^^, standing for the natural phase matching condition in the PCF 



for a monochromatic pump. The second one is e 



standing for the pump 



pulse broadening. So if we want to sum over all the possibilities for our fibre we have 
to look at the product of these two functions and to integrate it over Auo s and Aui. 
We plot on figure |U the graphical representation of both functions and their product. 
We can then carry out the formal calculation of (N s ) as explained in Appendix C to 
eventually find: 



(Ns) 




(21) 
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Figure 4. The hue color represents high values of the function. It is revealing to note 
the angle of the sine stripe which is controlled by the coefficients Af s ,Afi, N p — arising 
from the fibre geometry. This is slightly different from the angle of the exponential 
stripe which is always fixed at 45°. The difference in angle between the two functions 
defines the narrowness of the signal and idler bandwidth. 



If we now develop the Afc-function around the phase matched frequencies using 
first order Taylor power series but use the energy conservation to remove the Auj^ 
(respectively Auk s ) term, we get 



Awl 



Au k% 



2vrc 



\Ms-Ni\L 

2irc 
\K-Ns\L 



+ 2 



AT* - AC 



Aid- 



Auj T/ 



(22) 
(23) 



giving respectively the bandwidth of the signal and idler photons. It is worth noting 
two separated terms: the first one corresponding to the natural bandwidth from a 
monochromatic pump, while the second one is the pump bandwidth broadening. In 
the classical CW case, even a perfectly monochromatic pump leads to down-converted 
photon with a finite bandwidth. This process behaves as the inverse of the length of 
the fibre [22] ■ Here, we have to take it into account and add the fact that our pump is 
not monochromatic. The pump bandwidth contribution is linked to the slope 



for the signal and 



JV s -JVv 



Mi-Mv 



for the idler) of the phase matching curves in figure H3 



3.3. Experimental limitations 

In the previous treatment of the four-wave mixing process, we assumed that all the 
frequency components, and their following quantized operators, were all coherent within 
the pulse. Our calculation obviously applies to PCF pumped by Fourier transformed 
limited pulses, nevertheless it would be straightforward to decompose the pulse into 
several coherent subgroups, apply the treatment to and eventually sum up the number 
of photons coming from each of them. 

Another experimental limitation we did not highlight in the previous treatment 
is group velocity walk-off between the three involved photons. Although the pump is 
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close to the group velocity dispersion minimum, this is not the case for the signal and 
idler photons. After a certain distance, the pump pulse no longer overlaps with the pair 
pulse, thus meaning that we cannot sum up coherently all the operators. Using the 
propagation constant determined in section 13.11 we estimated the walk-off distance to 
be about lb cm (for our pulses). Therefore we have to apply our numerical calculation 
over this length and eventually to multiply accordingly to cover the actual length of the 
fibre. 

From (1221) and (j2HJ) , we should expect the bandwidths to remain constant beyond 
this length. Note however, in (|21|) . the number of pairs created always grows linearly 
with the length. As a consequence, the walk-off distance will not play a role in the 
number of created photons and monitoring the bandwidth fluctuation versus fibre length 
is the only way we have to experimentally estimate the figure associated to the walk-off. 



4. Experiment 

4-1. Setup 
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Figure 5. Optical layout. Laser, 708 nm Ti:Sa laser; P, prism; HWP, halfwave plate; 
PCF, photonic crystal fibre; M, protected silver mirror (R < 95%); DM, dichroic 
mirror (centered700 nm, T < 85%, R < 90%); Fl, 570 nm band-pass filter, bandwidth 
40 nm, T = 80%; F2, 880 nm band-pass filter, bandwidth 40 nm, T = 80%; APD, 
Silicon single photon detector (n s w 60% and i]i ss 33%). 



In order to estimate the brightness of our source we used the coincidence setup 
depicted in figure El where a mode-locked picosecond TkSapphire pump laser (Spectra 
Physics - Tsunami) emitting ~ 2ps pulses with a repetition rate of 80 Mhz is sent, 
through an optical isolator, onto a prism P to remove in-band light from the pump 
laser spontaneous emission. A pin hole is then used to improve the pump mode and 
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eventually several attenuators bring the power down so that up to 1 mW average power 
is launched into the fibre. Since the PCF is birefringent and supports two modes, a half 
wave-plate (HWP) is used to align the pump polarization along one axis thus preventing 
polarization scrambling and creating pairs with the same polarization as the pump beam. 
The output of the fibre is collimated using an aspheric lens, followed by a dichroic mirror 
centered at 700 nm to spread the incoming beam into two arms, one corresponding to 
the signal channel and the other to the idler, where band-pass filters Fl and F2 centered 
at 570 nm and 880 nm respectively (width 40 nm, T > 80%) are used to remove in-line 
pump and background light. Each photon of the pair is then launched into single mode 
fibres that are connected to two Silicon avalanche photodiodes (APD). The detected 
photons are counted in a dual-channel counter and the coincidences between the two 
APDs are analyzed using a time interval analysis system (TIA). 

4-2. Results 

Looking now at the coincidence test bench in order to determine the brightness of our 
source, we measured the number of single counts in both signal and idler channels, 
while we recorded the number of coincidences. This experimental protocol amounts to 
recording the percentage of detections in the "start channel" which have been stopped 
by detection in the "stop channel" in the following time interval and effectively provides 
a direct estimate of the "lumped" efficiency of the stop channel. 

4-2.1. Photon pair rate The central peak on figure El corresponds to signal and idler 
photons belonging to the same pulse. The small satellite peaks stand for uncorrelated 
events, i.e. signal and idler coming from subsequent pulses, whether they are actual 
pairs of photons or background photons. It is interesting to note the satellite peaks 
grow with pump power, whereas the central peak remains constant (the arrows point at 
its height). The central peak reflects the lumped collection efficiency which is relatively 
constant with increasing pump power (confirming correlated photon pair creation), while 
the satellite ones are linked to the probability of having a photon from a different pair or 
a background photon in a neighbouring pulse. In the absence of background the satellite 
peaks would thus be indicative of the rate of generation of pairs of pairs, or 4-photon 
events. It is these events that could be useful for multi-photon interference experiments 
and quantum information applications. However, we are here interested in producing 
no more than one pair of photons per pulse and we still have to take into account both 
the background and the multi-photon pair rate in the central peak as they increase the 
photon count probability. Thus we will introduce C raw as the raw coincidence rate in 
the central peak and we will calculate the accidental coincidence rate thanks to the 
satellites peaks C&. This allows us to write for C raw : 

Craw = ^sVs^iW + C b (24) 

where r is the actual photon pair rate; r] s and 77* are the APD quantum efficiencies at 
586 nm and 894 nm, while the net optical transmission and launch efficiencies into single 
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Figure 6. Time interval histogram showing the coincident photon detection peak and 
also a zoom on one of the accidental coincidence peak for different pump powers. The 
instrument displays the probability that a start pulse is stopped within a given time 
bin. Here the time-bin width is 156 ps. The time between two peaks reflects the pump 
laser repetition rate. However the width of the peaks is limited by the response time 
of the detectors which is typically hundreds of picoseconds (rather than the actual 
duration of the pump pulses). 

mode fibre of each arm are /i s and //j. We defined the single counting rates in the signal 
and idler APDs as: 

n=oo 

N r s aw = fx sr]s r + Msr n + B s (25) 

n=2 
n=oo 

Nl aw = mi r + £ ^ Vt r n + B t (26) 

n=2 

where B s , Bi are total background rates and fi s r] s r n , fiiT]ir n stand for the multiphoton 
pair contributions. Here, we separated the background and multiphoton contribution 
in ()25j) and (|26p. as we can already consider the latter being negligible for low average 
powers. For instance, as long as we remain under 0.1 photon pairs per pulse, the 
high order contribution should remain a tenth of the total. On the other hand, the 
background might not be negligible and is mostly due to Raman scattering in our 
case j2Hj- Assuming we can single out this background contribution, as shown in 
reference [31], we can now use the singles counting rates and the coincidence rates 
to estimate the actual rate of pairs r created inside the PCF using this equation: 

N r s aw = + B s 
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[S'raw ~ ^b) 

As previously stated, 5, and B s include the Raman background rates in the APDs 
and we have to single out these detections. The number of created photon pairs is 
proportional to the square of the peak intensity and the Raman scattering grows roughly 
linearly at these relatively low pump powers. Thus switching the laser from the pulsed 
regime to the CW regime reduces the photon pair rate to negligible levels while keeping 
the Raman and other sources of background constant for the same average power. We 
propose to validate the technique by recording the number of events in the two APDs 
versus the pump power for each case (i.e. pulsed and CW) and studying the evolution 
of the counting rate summarized in table ^ Here a linear tendency would be the proof 
of spontaneous Raman scattering, while pure quadratic behavior would be the proof 
of actual photon pairs. All the plotted data have been corrected to compensate the 
non-linearity of the Si-APD for high counting rate following the datasheet 



Table 1. Experimental measurement of the coincidence rate versus the pump power 
for 0.2 m of PCF. The dark count rate in both APD is - 400 Hz. 

Power N™ w [N° w ] N[ aw [Nf w ] C raw [C b ] 

960 fW 2070 kHz [5.0 kHz] 917 kHz [54 kHz] 234 kHz [16 kHz] 

660 uW 1031kHz [3.0 kHz] 460 kHz [38kHz] 115 kHz [3.8 kHz] 

490 uW 601kHz [2.0 kHz] 261kHz [25 kHz] 66 kHz [1.4kHz] 

340 iW 298 kHz [1.2 kHz] 135 kHz [18kHz] 33 kHz [0.1 kHz] 

200 uW 107 kHz [0.9 kHz] 52kHz [11 kHz] 12 kHz [0 kHz] 

0^W 0.4kHz [0.4kHz] 0.4kHz [0.4kHz] 0kHz [0 kHz] 



We first plot in figure the number of detections versus the pump power in CW 
regime. We clearly see the negligible level of Raman scattering at the signal wavelength 
compared to the idler wavelength, as expected as the stokes Raman scattering is 
often considered negligible compared to the anti-stokes one. In both cases, the fit is 
purely linear (quadratic term is negligible) and being the proof of spontaneous Raman 
scattering process and possible leakage of pump light. We can therefore use these figures 
as a good estimation of the Raman background rate for any given mean pump power. 

We then plot in figure |H1 the number of net detections versus the pump power. We 
define the net detections as being the number of detections in pulsed regime minus the 
number of detections for the same average power in CW regime. We clearly see that the 
number of net detection versus the pump power is purely quadratic, thus highlighting 
that we are dealing with actual photon pairs. We can then now identify B s and Bi as 
the CW counting rate and calculate in table |2] the actual number of photon generated in 
the fibre using (fTTj) . To perform the calculation, we don't have to correct the counting 
rate from the APD, include all kinds of losses for the twin photons (even 

APD non linearity). If we eventually compare the counting rate associated with the 
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Figure 7. Corrected raw counting rate as function of the pump power in both APD 
in CW regime for 0.2 m of fibre. 
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Figure 8. Corrected net counting rate as function of the pump power in both APD 
in pulsed regime for 0.2 m of fibre. The net counting rate stands for the Raw counting 
rate minus the CW counting rate. 
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actual Raman scattering to the counting rate associated with the actual photon pair, 
it confirms that the background in the signal channel is extremely low (< 0.3% of the 
total rate) while in the infra-red we estimate the background rate to remain relatively 
low (< 6%). 

Table 2. Experimental estimation of the photon pair rate versus the pump power for 
0.2 m of PCF. From left to right, we first recall the associated pump power and the 
theoretical photon pair rate using (|21|l . followed by the experimental photon pair rate 
using (|27J) . For information purpose, we added the mean number of photon per pulse 
( ft exp ), the coincidence to accidental contrast ( c ™™~ c >> ) and the lumped probability 
of detecting respectively the signal (fi s rj s ) and idler photons (jiigi). 



Power 


r th (s- 1 ) 


T'exp ) 


( Us )exp 


C/A 






960 /W 


~ 8.05 xlO 6 


8.46 xlO 6 


0.11 


15:1 


0.235 


0.106 


660 fW 


- 3.81 xlO 6 


4.08 xlO 6 


0.05 


36:1 


0.240 


0.109 


490 fW 


- 2.10xl0 6 


2.31 xlO 6 


0.03 


55:1 


0.244 


0.109 


340 iW 


- l.lOxlO 6 


1.14xl0 6 


0.015 


220:1 


0.244 


0.109 


200 iW 


- 0.35 xlO 6 


0.43 xlO 6 


0.006 




0.223 


0.109 



If we compare the photon pair rate to the theoretical rate, we have here good 
agreement. However, one has to bear in mind that the estimation of the mode area 
and the peak power could lead to large error in the expected pair rate rth- We found 
it useful to recall the lumped probability to detect the photons of the pairs. These 
important experimental parameters are useful for evaluating the possibility to perform 
multi-photon pairs experiment as discussed in section El Here, provided a detection 
efficiency r] s ps 0.60 and r\i rs 0.33 35], we found a global coupling efficiency from 
the PCF fibre to the single mode fibre of \i s ~ 0.42 and ~ 0.30. These figures 
include the fibre coupling, filters, dichroic mirror and lenses losses for signal and idler 
channels respectively. By correcting these results with the filters and dichroic mirror 
losses, we find a coupling efficiency into single mode fibre of respectively 0.58 and 0.44 
for the signal and idler photons. We also found useful to introduce the coincidence 
to accidental contrast. This figure of merit has been introduced in [36J to quantify 
the quality of any presented coincidence rate. Our results clearly illustrate that this 
measure can be misleading as the coincidence to accidentals ratio drops with increasing 
pump power. This is due primarily to the increasing number of overlaps between pairs 
of pairs, which we hope to exploit in future multiphoton interference experiments. 

4-2.2. Walk- off Taking into account the walk-off distance we identified in section l3~3*l 
we should expect the bandwidth to be constant versus the fibre length whenever the 
fibre is longer than 15 cm. For several fibre length, we recorded the pump bandwidth, 
measured the signal fluorescence bandwidth and compared it to (}2*2*j) . From figure |H1 
one can note the quite good matching between the measured and predicted bandwidth 
and the similar behavior of both bandwidths versus the fibre length. As predicted, the 
measured bandwidth remains constant whenever the fibre is longer than the walk-off 
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Figure 9. Comparison of the measured and predicted bandwidth of the emitted signal 
photons. 

distance which we estimate to be around 15 cm. More measurements around 15 cm would 
be needed to reduce the uncertainty about the walk-off experimental determination. 

5. Discussion 

This source has to be usable for future quantum interference experiments involving three 
or more photons created as two or more pairs [301 - W e dedicated the following section 
to estimate the potential of this source. To do so, we used realistic requirements for 
interference effects between separate pair-photons and attached experimental figures to 
a typical four-photon quantum experiment. 

5.1. Multi-photon pair weight in satellites peaks in figure^ 

We showed in the previous sections that the spontaneous Raman scattering was 
negligible compared to the photon pair rate. However this background may be equivalent 
in intensity to the probability of creating two pairs of photons. We can use the diagram 
in figure |H1 to state that our source exhibits a background low enough to allows multi- 
photon pairs experiments. The spontaneous Raman scattering rate grows linearly with 
pump power whereas the photon pair rate grows quadratically. Hence, in figure El we 
would therefore expect the satellite peaks to decrease in significance as we increase pump 
power as the ratio Raman/photon pair decreases. Experimentally, the satellite peaks in 
figure El increase clearly showing that the background coincidence rate at high counting 
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rates is dominated by random overlap of more than one pair of photons. The satellites 
peaks are thus linked to the multi-photon pair probability, they are indicative of the 
four photon coincidence rate we will get in future experiments and are variables we want 
to maximize. 

For these multi-photon quantum information experiments, we require a narrower 
bandwidth so that the coherence length is equivalent to the pulse length [3Uj. A quantum 
interference experiment involving four-fold coincidence between photons coming from 
two separated sources would require a filter of order 0.2 nm bandwidth in the green 
(0.4 nm in the IR). Such filters will transmit only 40% of in-band light thus halving 
our effective efficiencies and collect only ~ 1/25 of the available spectrum. However 
our present coincidence rates are limited by the detector saturation. Here the singles 
counting rates would be significantly reduced thus allowing an increase in pump power. 
Using our source with a pump power of ~ 6mW, the expected rate of photon pairs 
detected within this bandwidth is ~ 9 x 10 s -1 , which means a rate of four photon 
events ~ 100 s" 1 , two orders of magnitude higher than any previous experiment. It 
is also important to underline, that our Raman background rate in the APDs, will be 
reduced by a bandwidth factor of ~ 1/100 due to its broadband nature thus improving 
even more the performances of the source. 

6. Conclusion 

We have presented in this article both theoretical and experimental aspects of photon 
pair generation using four wave mixing in a photonic crystal fibre. We have clearly 
demonstrated a good agreement between the theory and experiment. Any realistic 
multi-photon quantum experiment requires a precise timing using ultra-short pulses, 
thus we have developed the first quantum model suitable for such a regime. We used it 
to suggest a numerical estimation of the photon pair rate from a photonic crystal fibre 
provided a good knowledge of the mode propagating in the fibre. We have validated the 
model by comparison with the experimental measurement of picosecond-pulsed photon 
pairs generated by four-wave mixing in a single-mode optical fibre. The source we used 
is polarized, bright, narrowband, single- mode and tunable by varying laser wavelength 
or fibre parameters. The wide separation of the generated pair wavelengths means that 
most of the background can be avoided. All these advantages make this new source of 
photon pairs more appropriated compared to conventional ones for quantum information 
processing applications. 
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Appendix A. Derivation of the interaction Hamiltonian in the pulsed case 

The interaction Hamiltonian in a volume within the fibre is: 

H mt = [ UalV (A.l) 



and considering the energy density associated with four fields coupled in a third order 
non-linear medium to be: 

U = e x (3) (Op: f „ uJi)EpE s Ei (A.2) 

In the four wave mixing-process, we will look at two pump frequency components 
{u p i,Up2} interacting one with each other. This new component of our interest is 

2 2 

a Gaussian defined as follow where G p (u) = e~ i± ^~ is the weight of one frequency 
component of the pump pulse: 

x G p (d p2 ) (e-^+^-^-^M + h.c.) e p2 {f) (A.3) 

E 2 

where -f 2 - is proportional to the density of energy per frequency units squared, and e p (r) 
is the transverse spatial variation with normalisation / \e p {r)\ 2 rdrdO = 1. We take the 
quantified fields for the signal (I = s) and idler (I = i) modes as: 



Ei = E (aje-^ - a ie ^) e t (f) (A.4) 

The four-wave mixing interaction only occurs if the energy conservation and phase- 
matching condition are satisfied: 

(2Q p + Dpi + Dpi) — d s — Di =0 (A. 5) 

(kpi + kp 2 ^ —kg — ki — 2-fPp ■ Hp = (A. 6) 

where {Dpi,D p2 } are the relative frequencies around the central pump fl p . Since we are 
in a collinear guided configuration, we can omit the vectorial terms by projecting on the 
z-axis. We finally get for eqs IA.5I and IA.61 

(2Q p + Dpi + Dpi) — D s — Di = (A. 7) 

(k pl + k p2 ) -k s -ki- 2 1 P P = (A.8) 

where ki >s>p are the wave-vectors (propagation constants) of the idler, signal and pump 
photons; P p is the peak pump power and 7 is the nonlinear coefficient of the fibre 

7=^- (A.9) 
\A ej} 

where n 2 ~ 2 x 10 -20 m 2 jW is the nonlinear refractive index of silica and A e jj is 
the effective area of the fibre mode. We then restrict the interaction Hamiltonian 
for one frequency component of the pump pulse to the phase matching situation 
when Ak pa 0. This approximation is justified since all the other combinations (i.e. 
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[2k p + k s — ki — 27P, . . . , 2k p + k s + fcj — 2 r yP]) of fields will lead to null terms after 
integration over the length of the fibre. 

r L/2 I 

where I = jj e p i(r^e p2 (r^e s (f)ei(r)drd9 is a normalised factor and Ak = k pl + k p2 — k s — 
ki — 2^P P the phase matching coefficient. It is convenient to define the gain parameter 
as: 



H int = SIh f L/2 \ G P M G p (u p2 ) f a l a \ e -^P+^+^)t e iAkz + h \ dz (A 1Q) 



E 



2 



By carrying out the integral over the length of the fibre, we get the following interaction 
Hamiltonian: 



+ h.c.) (A.12) 



It is important to note that the Hamiltonian 7i int only stands for a frequency component 
of the pulse. Therefore we will have to integrate later over the whole pump pulse 
frequencies. 

Appendix B. The Heinsenberg equation of motion 

In the standard Heisenberg representation the evolution of any operator A is defined by 
d i = -i\A,H M ] (B.1) 
In the case of A = a], equation IB. II gives 

^ = «*aj - £ G p (co pl ) G p M (a t e-^+^+^sinc j (B.2) 

We now assume that the gain S is small and that we can make the slowly varying 
envelope approximation (SVEA) ai(t) = ai(t)e~ iu ' lt where I = s,i leading to simplify 
eqlB~2linto: 

f = -^EG P M G P M (a^sinc (^Jj (B.3) 

where Auj = (2Q p + u p i + u> p2 ) — uj s — Ui. Since di(t) is slowly varying, we can assume 
that aj(t) ~ 5j(0) and for each Fourier components of the pulse, we can easily integrate 
over the time from to t —>■ 00 (we assumed that the system — nonlinear medium and 
pump laser — has been ready for a long time) leading to the 5-function: 

~a\(t -> 00) = al(0) - J2 G P (M G p (u p2 ) ^(0)8 (Auj) sine ) ( B - 4 ) 
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It is interesting to note here that the energy conservation lies in the ^-function of Au, as 
it would be in the monochromatic pump case. If we now sum up over the whole pump 
frequencies, we get: 

SIL 



JJ^2 G p ( u pi) G p (^2) ai(0)5(Au)sinc ( — — j du pl du p2 

It is possible to invert the integrals over {duj pl , du P 2} and the sum over fc, as they are 
independent variables. This leads us to integrate over a monochromatic field operator d\, 
for a given (monochromatic) mode in the idler. As a consequence the dirac function helps 
to compute the double integral since duj p \ has to move as the opposite of du p2 to keep 
satisfied the energy conservation. Given the basic energy conservation 2Q p = u SQ + u io , 
we are here performing the integration for a signal mode k s and an idler mode fcj 
satisfying the energy conservation 

(Q p + u pl ) + (ft p + u p2 ) = (uj so + Au k J + (u io + Au) ki ) (B.6) 

where we have introduced: 

• Au>k s — the frequency shift for the signal mode compared to the central one u So 

• Au ki — the frequency shift for the idler mode compared to the central one Ui 

For instance, (|B.6j) requires u P 2 to become u P 2 — Au^ (while u p i becomes uj p \ — Acj^J 
to satisfy the energy conservation. Eq IB. 51 becomes then: 



SIL 

ot(t-oo) =at(0)-i^x (B.7) 

rr _ Ki- A ^) 2 - 2 _ (^2-a^) 2 ^ (AkL\ 

2_^\\z 2 e 2 ai(0)o(Aui)sinc I du) p \du p 2 



k , 



\ 2 / 



The (5-function is used to reduce the number of integrals, by changing u P 2 into —uj p \. We 
can also introduce the new variable uj p \ — > cu[ — A " fcs 2 A ^ fc ' and assume that sine {^r) 
will not change significantly while we scan uj p \ over Auj p . We now find after integration 
over u[ from — oo to oo: 



12 



oJ(t - oo) = fit ( o) £ ^ e i a,(0)«nc I — J (B.8) 

In ()B.5|) . we were considering a monochromatic component of the pump pulse. Thus 
the theory still involved a Dirac delta function arising from the energy conservation. 
However, in (|B.8|) we have left the monochromatic pump regime and lost the Dirac 
function while integrating over the pump pulse. However it is interesting to note that 
we still have energy conservation in the broadband pump regime represented within the 
Gaussian function with width limited by the bandwidth of our pump pulse. In order 
to calculate the mean number of photons (N s ) over all the possible modes k s , we define 
the operator W s (t) to be given by a sum over the k s : 

W s (t) = ^=J2~ a s(t) (B.9) 
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so that the mean number of photon per pulse is given by: 

(N s ) =i (0| s <0|WtW s |0) s |0> i 
From Eqs. IB .81 and IB. 91 we found for (N s ) the following equation 
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a 



E E 



[a^+au^ „ 2 ( AkL\ 
4 sine I — : — I x 



V 2 / 



(B.10) 



(B.ll) 



sine 



f Ak'V 



<0| i (0|a i 4|0> s |0) i 



The commutation properties imply k s > = k s , k^ = ki and Ak' = Ak and allow us to 
simplify the formula 



w=(f9 ItIs 



h s Mi 



2 ( AkL\ hf^+^y ° 2 
sine e 2 

v 2 ; 



(B.12) 



In the limit L — > oo, the sum over k s ^ can be replaced by an integral in the usual way 

2 



2 /AfclA [^ ks +^] ° 2 
sine I — - — le 2 dk s dki (B.13) 



\ 2 j \ a 

As a final step, we can now convert the integrals over modes into integrals over 
frequencies [37] where v 9s and v gi are the group velocity of signal and idler wave in 
the fibre. 



(N s ) 



fSILy 



V 2 



■ 



7[ 



a 



sine 



/AkL\ \ A ^s+^} 2 * 2 dAu ks dAu, 
V 2 / % % 



(B.14) 



Appendix C. Estimation of the number of photons 

This integral can be numerically evaluated if we now develop the Afc-function around 
the phase matched frequencies using first order Taylor power series: 

A&; = — H p {dujpi + dujpz) + Af s du s + Afiduii (C.l) 
. Then with the help of energy conservation (—du p i — duj p 2 + 

(C.2) 



/ I, dn l 













where Mi 

duji + du s = 0) we find: 



^ _ "_JLL ] U j s _|_ Alii "JJLLdw, 



We can now write EqEH 

fSIL\ 2 



(N s ) 



V 2 



71 



X 



a 



(C.3) 



sine 



\N s -Np)L 



2c 



Au; fcs + 



(a/;-aql 



2c 



It is quite easy to integrate this product making the following change of variables 
(X = AuJks ^ AuJk i anc [ y = Au}ka AuJki ) . Then it becomes obvious that if we 
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perform the Y integral first and consider the sine function as a gate of height 1 and 
FWHM= (■ A/ " s ~ A/ '') L Ay = ir, the system is then easy to solve and we find: 
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